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ABSTRACT 

The 7 Cephei system is one of the most closely bound binary planet hosts known to date. The companion 
(7 Cep B) to the planet-hosting star (7 Cep A) should have truncated any protoplanetary disk around 7 Cep A, 
possibly limiting planet formation in the disk. We explore this problem by calculating the truncation radii of 
protoplanetary disk models around 7 Cep A to determine whether or not there is sufficient material remaining in 
the disk to form a planet. We vary the accretion rate and viscosity parameter of the disk models to cover a range 
of reasonable possibilities for the disks properties and determine that for accretion rates of > 10"^ Mq yr"' and 
low viscosity parameter, sufficient material in gas and solids exist for planet formation via core accretion to 
be possible. Disk instability is less favored, as this can only occur in the most massive disk model with an 
extremely high accretion rate. 

Subject headings: accretion disks — binaries: close — stars: individual (7 Cephei) — planetary systems — 
planetary systems: formation — planetary systems: protoplanetary disks 



1. INTRODUCTION 

Over the last decade, several stellar and substellar compan- 
ions of exoplanet host stars have been found, mostly in seeing- 
limited wide field i maging surveys (e.g., 'Mu grauer et alJ 
l2007bHRaghavan et al. 2006; Chauvin et al. 2006]). The pro- 
jected separations are a few tens to several thousa nd AU, with 
compa nion masses ^ 0.08- 1.1 Mq. Eggenberger & Udrvl 
(l2007h suggest that exoplanets may be less common in bina- 
ries closer than 120 AU. This apparent lack of close compan- 
ions to exoplanet host stars may indicate that planet formation 
is hampered by the gravitational influence of a close mas- 
sive companion. However, further investigations are needed 
to confirm this result. 

One of the closest planet host binaries presently known 
is 7 Cep. The planet orbits the prim ary 7CepA on a 906 
day orbit and has msin(0 r^ \.l Mjup ( Campbell et al.llT988l : 
iHatzes et al.ll2003h . The stellar companion orbi ts at a semi- 
major axis of 20 A U with an eccentricity of 0.4 (iHatzes et al.l 
l2003t lTbrresll200'7n . Given that massive companions can dis- 
rupt the protoplanetary disks in which planets form, how fea- 
sible is the in situ formation of a planet around 7 Cep A? 

Haghighipour (2006) considered the dynamical stability of 
7 Cep in order to put limits o n possibility of the e xistence 
of another planet in the system. iThebault et al.l ( |2004|) consid- 
ered the problem of planet formation via core accretion in this 
system using N-body simulations, assuming a density profile 
consistent with the Minimum Mass Solar Nebula (MMSN) 
model of iHayashil (1 198 ID . which is ste eper than that pro- 
duced by viscous accretion disk models. iPaardekooper et aP 
(I2OO8I) revisited the problem including gas drag and deter- 
mined that giant planet formation by core accretion is feasible 
in 7 Cep, although they did not address the evolution of the 
gas in the disk or the truncation radius of the disk in detail. 
iKlev & NelsonI (12008;) assumed a specific disk model and ex- 
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amined the fate of planet cores inserted into the disk. Rather 
than do a detailed hydrodynamic simulation of a planet em- 
bedded in the disk, our objectives are to model the protoplan- 
etary disk in 7 Cep as a viscous accretion disk and to explore 
which disk parameters allow planet formation to occur given 
that the disk is truncated by the stellar compa nion. 

This analysis is similar to that done in Jang -Condelll (l2007h 
(henceforth Paper 1) for the extremely close triple sys- 
tem HD 188753. Paper 1 concluded that HD 188753 was 
unlikely to support a disk sufficiently massive to support 
planet formation. Indeed, the initial clai m of a Jupiter- 
mass plane t in HD 188753 (JK onacki 2 005|) has since been 
refuted (Eggenb erger et all 1200 7). This does not rule out 
the possibility that the planet could form around a single 
star or in a wide binary and then undergo dynamical evo- 
lution, such as through cl ose encounters with another star 
(iPortegies Zwart & McMil lan 2005. lPfahll[2005h . but this is 
outside the scope of this paper. 

2. MODEL DESCRIPTION 

We adopt orb i tal pa rameters for the 7 Cep system from 
iNeuhauser et al.l (l2007h . as follows: primary mass I.4OM0, 
secondary mass 0.409 M©, eccentricity 0.41, and semi-major 
axis 20.18 AU. We ignore the orbit of the planet, since we 
are interested in pre-planetary conditions of the disk around 
the 1 .4 Mq primary. We assume that the stars have not un- 
dergone significant mass loss or orbital evolution since their 
formation, so we can model the properties for 7 CepA based 
on a pre-main sequence stellar model for a 1 .4 Mq star. Since 
the typical age of a T Tauri star is 1 Myr, we assume this age 
for our model. 

2.1. Disk Model 

The calculation f or the disk models is de s cribed in de - 
tail in Paper 1 and iJang-Condell & Sasselovl (l2003l |2004|) . 
We assume an a-disk model, where the viscosity v is given 
hy V = acji where Cj is the sound speed, h is the thermal 
scale height of the disk, and a is a dimensionless parameter 
{Shakura & Sunvaev 1973; Pringle 1981). The disk temper- 
ature is set by stellar irradiation at the surface and viscous 
heating at the midplane. The radial and vertical density and 
temperature structure of the disk are calculated iteratively for 
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self-consistency. We adopt effective temperature T* = 4500 
K, and radius /?* = 3.0 R©, correspondi ng to a M^, = 1.4M ©, 
1 Myr old star with metallicity Z = 0.02 dSiess et al.ll2000h . 

The two remaining free parameters for our disk models 
are the mass accretion rate onto the star M, and the vis- 
cosity parameter a. The exact values for these parameters 
are unknown for 7 Cep, and it is likely that these values 
evolved over time, so we explore a range of values for both 
M and a to determine which, if any, set of parameters al- 
lows for planet formation to occur As in Paper 1, we cal- 
culate a grid of disk models, with a G {0.001,0.01,0.1} and 
M e {lO-^lO-^lO-^lO-^lO-^lO-'*}M0yr-l. These pa- 
rameters are roughly consistent with observations of T Tauri 
stars (e.g. Gullbring et al. 19981: iHartmann etall 11998) . in- 
cluding the extreme ly high and transient accretion rates of 
FUOri phenomena (ICalvet et alJl2000t IHartmann & KenyonI 
Il996l) . In practice, the models are calculated out to 256 
AU, but we consider only the material interior to the trun- 
cation radius to be available for planet formation. We re- 
fer to a given disk model by the coordinate pair (a,M), so 
that Model (0.01, 10"^) refers to the run with a = 0.01 and 
M=10-^MQyr-i. 

2.2. Disk Truncation 

The truncation radii of gaseo us disks depend on the viscos- 
ity and temperature of the gas ("Artymowicz & Lubowll 19941 
henceforth AL), as opposed to planetesimal disks, whose 
truncation radii can be calculated from last stable orbits of 
test particles ( Pichardo et al.ll2005h . The truncation radius of 
each disk model is calculated following AL, as in Paper 1. 
In AL, the truncation radius of a circumstellar disk in a close 
binary is where resonant and viscous torques balance. This 
depends on the mass ratio of the binary (fi), the semimajor 
axis of the orbit (a), the eccentricity of the orbit (e), and the 
Reynolds number of the disk (Re). For 7 Cep, /i, a, and e have 
all been determined observationally. The remaining parame- 
ter. Re, depends on the structure of the disk, which has long 
since dissipated. 

The Reynolds number is defined as Re = rv^jv where r is 

distance from the star, v^ = yjGM^jr is the orbital velocity, 
and V = acji is the viscosity of the disk. Since q/v^ = hjr 

and Cs = y/kT/m, 
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Setting e = 0.41 for 7 Cep, we read off truncation radii in 
units of semi-major axis of the circumprimary disk versus Re 
from Figs. 5 and 6 in AL, for fi = 0.3 and 0.1, respectively. 
For 7 Cep, /i = 0.22, so we interpolate in /i to find the final 
truncation radius versus Re relation. This relation is plotted 
as long-dashed black line in Fig. |2] for a semi-major axis of 
20.18 AU. 

We assume that the disks are dynamically truncated and that 
irradiation from the stellar companion is negligible compared 
to heating from viscous accretion and the central star. This 
irradiation would most likely further decrease the likelihood 
of planet formation since it would provide an additional heat 
source at the outer edge of the disk, inhibiting planet forma- 
tion by either core accretion or disk instability. In the absence 
of additional accretion of material past the companion's or- 
bit onto the disk, the disk should be viscously spreading both 
inwards and outwards. Thus, the calculated truncated disk 
masses should be considered upper limits. 




Fig. 1 . — Enclosed disk mass versus radius for dislc models for 7 Cep, in 
units of Mj. The accretion rate is indicated by color and symbol type: orange 
asterisks for 10^, red circles for 10"', green triangles for 10"*, blue squares 
for 10"', cyan stars for 10"^, and magenta crosses for 10"' Mq yr"' . Models 
with Q of 0.001, 0.01 and 0.01 are indicated by solid, dotted and dashed 
lines, respectively. The locations of the points mark the truncation radius and 
maximum disk mass for each disk model. 

3. RESULTS 

In Fig. [T] we plot the mass profiles of the disk models. The 
line colors and types (solid/dotted/dashed) indicate M and a 
parameter for each disk model, respectively (see legend for 
details). Disk mass increases with increasing M and decreas- 
ing a. The truncation radius for each disk model is indicated 
by a symbol on the line, which also indicates the maximum 
disk mass for each model. 

Fig. |2] shows Re versus radius for each disk model, coded 
the same as in Fig. [T] The Reynolds number in each disk 
model depends on the input parameters, but stays fairly flat 
with radius. The black long-dashed line shows the truncation 
radius versus Re relation calculated following AL. From the 
intersection of the long-dashed line with each model profile, 
we determine a unique truncation radius for each disk model, 
marked by symbols. 

The truncation radii are in the range of 4-7 AU, roughly 
consistent with the 4-5 AU truncation radius assumed by 
iThebault et al.l (l2004t) . albeit a bit larger. This is expected be- 
cause viscous torques of a gaseous disk allow it to extend far- 
ther than a particle-only disk. If we consider 1 .6 My, the mass 
of the planet 7 Cep Ab, to be the minimum mass necessary for 
in situ planet formation, a majority of the disk models satisfy 
this criterion. 

4. DISCUSSION: CORE ACCRETION VERSUS DISK INSTABILITY 

Having determined that a truncated disk around 7 Cep A 
can easily contain enough material to form one or several 
Jupiter mass planets, we now turn our attention to whether 
planet formation may take place via core accretion or gravita- 
tional instability. 

Giant planet formation by core accretion requires sufficient 
mass of solid material to coagulate into a dense core which 
can then accrete gas. We calculate the mass of dust or solid 
materials as in Pape r I, adopting the dust composition from 
Pollac k et al.l (Il994h of olivines, orthopyroxene, iron, water, 
troilite, refractory organics and volatile organics. At a given 
radius in the disk, we assume that each species is completely 
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Fig. 2. — Reynolds number versus radius for each of the disk models, 
marked the same as in Fig.[T] The black long-dashed line shows the truncation 
radius versus Reynolds number relation for the 7 Cep binary. The truncation 
radius for each model is gX the intersectiop with this line 

condensed or vaporized, depending on whether the tempera- 
ture is lower or higher, respectively, than the sublimation tem- 
perature for that species. 

In Fig. [3] we plot the mass of solids or dust particles in the 
disk as function of radius. Here, we only plot those disks that 
contain more than 1 .6 Mj of material within their truncation 
radii. We use the same coding of colors, lines and symbols 
to represent the different disks as in Fig. [U scaling the sym- 
bols sizes relative to the disk masses. The amount of solids 
does not simply scale with the mass of the disk because hotter 
disks sublimate their dust. While higher accretion rates yield 
more massive disks, they also yield higher temperatures. On 
the other hand, higher a values yield less massive disks and 
higher temperatures. If we require a minimum of 20 Earth 
masses (M®) of solids to form a giant planet core, we find 
that the Models (0.001, 10"^), (0.01, lO'^), (0.001, 10"^), and 
(0.001 , 10"^) all have sufficient solids to form giant planets by 
core accretion. 

The metric for planet formation by gravitational instability 
is the Toomre Q parameter: 



Q- 



ttGE 



(2) 



In order for planet formation to proceed, Q < 1 is required. 
In Fig. |4]we plot the local value of Q versus radius for our 
disk models. Only six of these have Q < 10 and Model 
(0.001, 10"^), the most massive disk, has Q < 1. As discussed 
previously, 10~^MQyr"' is an extreme accretion rate, seen 
only episodically in FU Ori stars. Moreover, Q < I only in 
very the outermost part of the disk. Therefore, we find it un- 
likely that 7 Cep Ab formed through gravitational instability. 
These results may be further generalized to suggest that 
planet formation in any close binary is more likely to occur 
through core accretion than disk instability. This is because 
disk instability happens only in the most massive disks. While 
the amount of solids is also dependent to some extent on the 
mass of the disk, the minimum disk mass for core accretion is 
still below that required for disk instability. 

5. CONCLUSIONS AND FUTURE WORK 

We assessed the feasibility of in situ planet formation in 
7 Cep by examining the properties of a protoplanetary disk 
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Fig. 3. — Mass of solids for each disk model. The position of the point 
indicates the truncation radius and maximum enclosed solid mass. The sizes 
of the points indicate the relative total disk mass. See Fig.[T]for key. 

around the primary star, given the current orbital parameters 
of the binary star. We examined a range of accretion rates 
and viscosity parameters and determined the truncation ra- 
dius for each disk model. We find that 7 Cep A can host a 
truncated disk of sufficient mass to form a giant planet with 
reasonable accretion rates and viscosity parameters, so in situ 
planet formation is possible. There are sufficient solids in the 
truncated disk for core accretion to occur for accretion rates 
higher than lO~^M0yr"' and low values of the viscosity pa- 
rameter, a. This is a relatively high accretion rate, which may 
indicate that giant planet formation must take place very early 
on, within 10^-10^ years, since it appea rs that accretion rates 
of T Tauri stars decrease with age (e.g. I SiciUa-Aguilar et al] 
I2OO6I) . Giant planet formation by disk instability is unlikely 
to have occured in 7 Cep. Disk instability also requires an 
extremely high accretion rate, ^ 10~^MQyr"', a rate that is 
typical of a transient FU Ori phenomenon. On the other hand, 
this may mean that FU Ori outbursts and giant planet forma- 
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Fig. 4. — Toomre Q parameter for each disk model, as indicated by the 
labels. The point on each line indicates the truncation radius. Models not 
shown on this plot have Q > 10 interior to their truncation radii. 
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tion by disk instability are correlated. 

We have omitted effects such as shock-heating or triggered 
plane t formation by the binary. Whether these effects inhibit 
(e.g.. lNelsonll2000l) or enhance (e.g.,IBoss 2006) planet for- 
mation depends on whether cooling times are long or short, 
respectively (see also Mayer et al. 2005) . Low mass disks that 
are not self-gravitating are not subject to s trong shocks, so ou r 
results for core accretion are unaffected dMaver et al.ll2005b . 
This strengthens our argument that core accretion is the fa- 
vored mechanism for giant planet formation in close binaries. 

The analysis presented here should hold also true for wider 
binary systems. That is, giant planet formation is possi- 
ble around any star where the stellar companion has a mass 
ratio // < 0.2 and has an orbit wider than 7 Cep B, with 
disk instability becoming increasingly feasible in wider sys- 
tems. The handful of triple systems with exoplanets that have 
been discovered to date are all hierarchical, where the exo- 
planet host star and a close stellar pair revolve around a com- 
mon barycenter. In these cases, the close stellar pair can be 
treated dynamically as a single object. The closest planet 
host triple system with confirmed ex oplanets presently know n 
is HD 65216 Ah-BC, at '-250AU (iMugraueret all i2007bl) . 
Thus, we conclude that there is no barrier to planet formation 
in th e known planet-hosti ng triple systems. 

In iJang-Condelll (l2007l) . the same analysis was carried out 
for a hypothesized planet in HD 188753, a multiple system 
which at first glance is only slightly closer than 7 Cep: a semi- 
major axis of 12.3 AU and eccentricity of 0.5. In situ planet 
formation in HD 188753 was ruled out, whereas it is deemed 



feasible for 7 Cep. Somewhere between the parameters of 
these two systems lies the transition between possiblity and 
impossibility of planet formation in close binaries. We will 
explore this parameter space and put limits on the closeness 
of planet host binaries in a future paper. 

Another interesting case is a system where the massive 
companion is a white dwarf. Known white dwarf planet 
hosts include G186B with a projected separation of only 
20 AU; and HD 2 7442 with a subgiant primary and white 
dwarf secondary (Mugra uer et alj |2007a). Stellar evolution 
and mass loss generally widens the orbits of binary stars, so it 
is likely that these systems were originally much closer (e.g. 
iDebes & Sigurdssonll2002l) . In order to determine whether in 
situ planet formation could have taken place in these sytems, 
the original orbital configuration must first be determined and 
the disk truncation radii based on that. We will address this 
issue in a future paper 
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